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Statement  of  the  problem 

The  objectives  of  this  research  were  to  obtain  spectra  from  impact  collisions  in 
the  visible  and  infrared  spectra  on  military  targets  that  would  collide.  The  Computed 
Tomographic  Imaging  Spectrometers  (CTIS)  can  hyperspectrally  characterize  transient 
events  in  an  imaging  format.  Temporal  changes  that  cause  noise  in  classical  scanning 
imaging  spectrometers  are  captured  by  the  CTIS  as  signal.  Tactically,  Battle  Damage 
Assessment  (BDA)  benefits  the  most  from  CTIS’s  ability  to  characterize  spatial  position 
and  spectrum  of  transient  events  which  were  the  objectives  of  these  measurements  made 
at  White  Sands  Missile  Range  (WSMR). 

The  University  of  Arizona  developed  a  state-of-the-art  mid-wave  infrared 
(MWIR)  and  visible  imaging  spectrometer  to  accomplish  these  tasks.  These  high-speed, 
non-scanning  Computed  Tomography  Imaging  Spectrometers  were  be  retrofitted  to  an 
existing  tracking  mounts  belonging  to  the  Army.  Tests  were  conducted  in  WSMR  to 
demonstrate  this  unique  capability. 


Summary 

The  driving  concept  behind  the  development  of  the  computed  tomographic 
imaging  spectrometer  (CTIS)  is  the  reconstruction  of  a  3D  object  from  2D  projections. 
The  objective  here  is  to  record  2D  projections  of  the  x,  y,  X  objeet  cube.  If  sufficient 
projections  are  recorded,  then  the  original  object  cube  can  be  reconstructed.  These 
projections  are  shown  schematically  in  Fig.  1.  A  two  dimensional  array  records  these 
many  “shadows;”  and  reconstruction  of  the  (x,  y,  X.)  object  is  accomplished  via  inversion 
math  methods. 


Fig.  1 .  Projections  of  the  object  cube.  Spectral  information  along  the 
vertical  axis  of  the  object  cube  is  projected  along  the  radial  coordinate  of 
the  focal  plane  and  multiplexed  with  spatial  information. 
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Description  of  the  Instrument 

The  spectrometer  consists  of  three  optical-element  groups;  an  objective  lens,  a 
collimator  lens,  and  a  re-imaging  lens  both  the  visible  and  infrared  instruments  are 
similar  in  optical  layout. 


Instrument  Layout; 

Objective  Lens 


Disperser 


FPA 


Collimating  Re-imaging 
lens  lens 


Fig.  2.  Optical  schematic  of  the  CTIS 


The  use  of  AR  coated  lenses  for  the  3  -  5  |j.m  band  throughout  the  optical  train  greatly 
enhances  the  performance  of  the  instrument.  The  Computer  Generated  Hologram  (CGH) 
disperser  is  located  in  collimated  space  between  the  collimator  lens  and  the  re-imaging 
lens.  A  250  mm  focal  length  objective  lens  focuses  the  target  onto  a  5mm  square  field 
stop,  thus  a  20  milliradian  field  of  view.  The  field  stop  is  appropriately  minified  and 
dispersed  on  to  the  InSb  FPA  by  the  combination  of  a  200  mm  collimating  lens,  GaAs 
binary  diffraction  grating  and  50  mm  re-imaging  lens.  The  FPA  in  the  prototype  system 
shov/n  in  Fig  3.  is  a  160  x  120  InSb  array  with  50  |im  pitch  and  30%  fill  factor  operating 
at  30  frames  per  second  and  5  msec,  integration  time.  The  Stirling  cycle  cooler  built  in  to 
the  camera  enhances  portability  of  the  system.  The  5  mm  square  field  stop  maps  to  an  11 
X  11  pixel  area  on  the  FPA.  Only  the  four  first  orders  are  used  in  reconstruction. 
Datacubes  up  to  1 1x11x8  (\,y,X)  have  been  reconstructed,  requiring  approximately  30ms 
per  reconstruction  on  a  portable  computer  (PII-333).  The  visible  instrument  uses  a  zoom 
lens  objective  with  a  5  mm  square  field  stop  have  a  20  milliradian  field  of  view.  The 
collimator  lens  (efl=  210  mm)  and  re-imaging  lens  (efl=28  mm)  makes  an  80  x  80  pixel 
area  with  32  spectral  channels  from  420  to  740  nm. 


Fig  3.  Prototype  Visible  CTIS  (showing  optical  train)  and  a  MWIR  CTIS 
inside  the  white  enclosure. 

While  the  ER  lenses  were  purchased  from  vendors,  the  two-dimensional  binary 
disperser,  which  is  actually  a  diffractive  phase  grating,  shown  in  Figure  4  was  custom 
fabricated  at  the  Optical  Sciences  Center  in  GaAs.  The  phase  grating  has  a  period  of  90 
|im  in  both  the  x  and  y  directions  was  chosen  in  combination  with  the  field  stop  size,  the 
collimating  lens  focal  length  and  the  re-imaging  lens  focal  length  to  place  the  zero  and 
first  diffraction  orders  with  in  the  focal  plane.  The  physical  depth  of  the  grating  was 
initially  chosen  to  be  0.7  |j.m.  The  depth  was  subsequently  fine  tuned  by  trial  and  error  to 
1 . 1  |J.m  to  give  reasonable  irradiance  uniformity  across  the  focal  plane. 


842.77 


Fig.  4.  Wyko  interferometer  surface  profile  of  the  2-D  binary  phase  grating  fabricated 
using  Reactive  Ion  Etching  at  the  Optical  Sciences  Center  by  Michael  Descour  and 
Daniel  Simon. 

CALIBRATION 

The  visible  system  is  calibrated  using  a  fiber  optics  and  a  Jarrel  Ash  grating 
spectrometer.  The  resolution  of  the  instrument  is  limited  by  the  calibration  process  to  be 
10  nm.  Each  resolution  element  is  calibrated  spectrally  over  the  visible  region  (80  x80 
x32). 

The  MWIR  CTIS  calibration  process  accomplished  using  a  973  K  blackbody  and 
a  grating  monochromator.  The  output  port  of  a  973K  blackbody  is  coupled  directly  into  a 
computer-controlled  monochromator  with  a  4-|im  blazed  grating.  The  monochromator 
allows  selection  of  spectral  bands  as  narrow  as  0. 1  |J.m.  The  output  of  the  monochromator 
must  be  imaged  into  the  field  stop  (via  a  fold  mirror)  from  a  distance  using  the 
spectrometer’s  objective  lens.  The  exit  slit  of  the  monochromator  is  stopped  down  to  a 
square  with  an  opaque  aperture. 

The  fold  mirror  is  used  to  position  the  calibration  spot  in  the  field  stop  so  that  it 
exactly  fills  a  single  (50-|im)  pixel  on  the  CCD.  The  monochromator  is  then  scanned 
through  wavelengths  and  a  calibration  image  is  acquired  at  each  spectral  band.  Since  the 
instrument  is  essentially  shift  invariant  over  each  wavelength,  the  calibration  images  may 
be  software-shifted  to  fill  the  field  stop.  There  is  currently  no  allowance  for  absorption 
over  the  path  difference  external  to  the  objective,  but  absorption  over  such  a  short  path 
should  not  be  significant  outside  of  the  CO2  absorption  band  from  4.2  to  4.4  )j,m.  There 
is  also  no  normalization  of  the  calibration  images  by  the  irradiance  of  the  incident  light. 
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Field  Tests 


The  field  tests  were  conducted  at  WSMR  on  KTM  tracking  mount.  A  photograph  of  the 
instrument  mounted  on  the  tracker  is  shown  in  Figure  5.  The  field  of  view  of  the 
instruments  were  within  the  tracking  accuracy  of  the  KTM  tracker,  however  there  was 
some  concern  about  keeping  the  target  within  the  instruments’  view. 


Figure  5.  CTIS  system  on  the  KTM  tracker  at  WSMR. 

The  calibration  data  on  a  819  K  blackbody  in  the  MWIR  is  shown  in  Figure  6.  The 
Blackbody  was  at  a  range  of  450  ft.  and  the  spectral  radiance  shows  a  dip  in  the  4.28  pm 
region  due  to  the  CO2  absorption.  Since  the  bands  are  wider  than  the  CO2  line  it  doesn’t 
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drop  to  zero  as  expected.  On  the  right  side  of  Figure  6,  the  spectral  images  are  shown  at 
different  wavelengths.  Again  the  4.25  |a.m  image  has  the  carbon  dioxide  feature. 


3.5  4.0 

Wavelength  (Microns) 


■ 

J.Obpm 

■ 

Ub 

■ 

3.4b 

■ 

3.6b 

■ 

3.8S 

m 

4.0b 

425 

■ 

4.45 

Figure  6  a.  MWIR  CTIS  looking  at  a  819 
K  blackbody  450  feet.  Dip  at  4.25  is  due  to 
CO2. 


Figure  6  b.  MWIR  spectral  images  of  a 
819  K  blackbody. 


Results 

We  have  demonstrated  a  high-speed  snapshot  imaging  spectrometer  prototype 
operating  in  the  MWIR  as  well  as  the  visible  in  a  field  deployable  portable  setup.  The 
instrument  is  currently  undergoing  modifications  for  increased  spatial-spectral  resolution 
and  for  field  portability.  Future  work  includes  design  of  a  more  efficient  disperser  and 
installation  of  a  large  format,  variable  integration  time  camera.  The  resulting  data  was 
delivered  to  WSMR  personnel. 
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ABSTRACT 

We  report  results  of  experimentation  with  a  new,  high-resolution  MWIR  non-scanning,  snapshot 
imaging  spectrometer  capable  of  simultaneously  recording  spatial  and  spectral  data  from  a  rapidly 
varying  target  scene.  The  instrument  is  based  on  computed  tomography  concepts  and  operates  in  a  mid¬ 
wave  infrared  band  of  3.0  to  5.0  |a,m.  High  speed  spectral  imaging  was  demonstrated  by  collecting 
spectro-spatial  snapshots  of  filtered  blackbodies,  combustion  products  and  targets  of  opportunity  in  the 
lab.  Raw  images  were  recorded  using  a  512  x  512  InSb  focal  plane  array  in  snapshot  mode. 

In  addition,  a  complete  Stokes  imaging  polarimeter  for  the  3-5  |im  band  has  been  constructed 
and  tested.  The  system  features  a  rotatable  form-birefringent  retarder,  wire  grid  polarizer  and  a  512  x 
512  InSb  array.  Design  considerations  are  discussed  and  results  of  experiments  with  man-made  and 
natural  targets  are  presented. 


MWIR  Snapshot  Imaging  Spectrometer: 
Calibration  and  Imaging  Experiments. 


INTRODUCTION 

Collection  of  spatial  and  spectral  information  has  both  military  and  commercial  applications, 
where  not  only  the  location  of  the  event  is  important  to  the  observer  but  also  the  spectral  content  for 
identification.  Up  to  now,  simultaneous  imaging  and  spectrometry  were  impossible.  Spectral  images 
were  obtained  by  scanning  in  time  or  space.  Recent  research  developments  have  produced  an  imaging 


spectrometer  which  simultaneously  collects  both  spectral  and  2D  spatial  information  about  a  scene.'  To 
demonstrate  our  concept,  it  helps  to  think  of  a  hyperspectral  object  cube  as  shown  in  Fig.  1. 


The  object  cube  represents  the  x,  y,  X.  information  of  a  scene.  Any  slice  through  the  cube  parallel 
to  the  x-y  plane  represents  a  monochromatic  image.  Shown  in  Fig.  1  are  different  conventional 
approaches  to  acquiring  the  data  within  the  object  cube.  Note  that  they  all  require  SCANNING  to  fill  the 
three  dimensional  space.  The  Computed  Tomographic  Imaging  Spectrometer  (CTIS)  obtains  information 
about  the  entire  object  cube  each  frame  time  as  shown  in  Fig.  2.  This  allows  reconstruction  of  spatial  and 
spectral  information  for  the  scene  at  each  frame  time,  which  is  essential  for  scenes  that  change  in 
space/wavelength  with  time.  The  CTIS  approach  to  hyperspectral  imaging  was  facilitated  by  the 
development  of  large  two-dimensional  arrays  and  fast  computing  capabilities.  In  order  to  acquire  the 
spectra  of  the  scene  on  a  broad-band  detector  array,  an  imaging  technique,  based  on  computer  axial 
tomography  (CAT),  was  adopted.  This  combination  of  technologies  allowed  the  development  of  a  no 
moving  parts,  optically  simple  instrument  called  CTIS.^ 


OVERVIEW  OF  CTIS  INSTRUMENT 

The  driving  concept  behind  the  development  of  the  CTIS  is  the  reconstruction  of  a  3D  object 
from  2D  projections.  The  objective  here  is  to  record  2D  projections  of  the  x,  y,  X  object  cube.  If 
sufficient  projections  are  recorded,  then  the  original  object  cube  can  be  reconstructed.  These  projections 
are  shown  schematically  in  Fig.  3.  A  two  dimensional  array  records  these  many  “shadows”  and 
reconstruction  of  the  (x,  y,  X)  object  is  accomplished  via  inversion  methods. 


Fig.  3:  Projections  of  the  object  cube  onto  the  focal  plane.  Spectral 
information  along  the  vertical  axis  of  the  object  cube  is  projected  along 
the  radial  coordinate  of  the  focal  plane  and  multiplexed  with  spatial 
information. 


The  spectrometer  consists  of  three  optical-element  groups:  an  objective  lens,  a  collimator  lens,  and  a  re¬ 
imaging  lens. 


Objective  Lens 


Instrument  Layout 

Disperser 


Fig.  4:  Optical  schematic  of  the  CTIS 


The  use  of  AR  coated  lenses  for  the  3  -  5  |i.m  band  in  the  optical  train  greatly  enhances  the 
performance  of  the  instrument.  The  Computer  Generated  Hologram  (CGH)  disperser^  is  located  in 
collimated  space  between  the  collimator  lens  and  the  re-imaging  lens.  The  camera  in  the  second 
generation  system  shown  in  Fig  5.  uses  a  512  x  512  pixel  InSb  array  with  25  |a,m  pitch  operating  at  60 
frames  per  second  and  12  msec,  snapshot  integration  time.  The  5mm  square  field  stop  maps  to  a  84  x  84 
pixel  area  on  the  Focal  Plane  Array  (FPA). 


Fig  5:  Second  generation  MWIR  CTIS.  From  left  to  right,  a  299  mm 
focal  length  objective  lens  focuses  the  target  onto  a  5mm  square  field 
stop.  The  field  stop  is  appropriately  minified  and  dispersed  on  to  the 
InSb  FPA  by  the  combination  of  a  250  mm  collimating  lens,  GaAs 
binary  diffraction  grating  and  100  mm  re-imaging  lens.  The  re-imaging 
lens  is  attached  to  the  liquid  nitrogen  camera  dewar. 


While  the  IR  lenses  were  purchased  from  vendors,  the  two-dimensional  binary  disperser,  which  is 
actually  a  diffractive  phase  grating,  was  custom  fabricated  at  the  Optical  Sciences  Center  in  GaAs.  The 
phase  grating  shown  in  Fig.  6,  which  has  a  period  of  90  pm  in  both  the  x  and  y  directions,  was  chosen  in 
combination  with  the  field  stop  size,  the  collimating  lens  focal  length  and  the  re-imaging  lens  focal  length 
to  place  the  zero  and  first  diffraction  orders  with  in  the  FPA  periphery.  The  physical  depth  of  the  grating 
was  initially  chosen  to  be  0.7  pm  based  on  a  model  derived  from  Goodman’s  treatment  of  a  sinusoidal 
phase  grating."*  The  depth  was  subsequently  fine  tuned  by  trial  and  error  to  1 . 1  pm  to  give  accceptable 
irradiance  uniformity  across  the  FPA. 
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Fig.  6:  Wyko  interferometer  surface  profile  of  the  2-D  binary  phase 
grating  fabricated  using  Reactive  Ion  Etching  at  the  Optical  Sciences 
Center  by  Michael  Descour  and  Daniel  Simon. 


CALIBRATION 


The  MWIR  CTIS  calibration  process  was  similar  to  the  visible  CTIS  calibration  process.  The 
output  port  of  a  commercial  blackbody  was  coupled  directly  into  a  computer-controlled  monochromator 
with  a  4-pm  blazed  grating.  The  monochromator  allowed  selection  of  spectral  bands  as  narrow  as  0. 1 
pm.  The  output  of  the  monochromator  was  collected  by  a  36  mm  F/1.6  compound  lens  and  launched 
down  a  0.5  meter  As-Se-Te  fiber.  The  fiber  250  pm  diameter  fiber  output  face  was  minified  into  the 
■nstrument  field  stop  with  a  25  mm  Ge  lens. 

The  Ge  minifying  lens  was  used  to  position  the  calibration  spot  in  the  field  stop  so  that  nine  (25- 
pm)  pixels  were  filled  on  the  FPA.  The  monochromator  was  then  scanned  through  21  wavelengths  and  a 
calibration  image  was  acquired  at  each  spectral  band.  Since  the  instrument  is  essentially  shift  invariant 
over  each  wavelength,  the  calibration  images  were  software-shifted  to  fill  the  field  stop.  This  technique 
used  to  calibrate  the  MWIR  CTIS  with  28  x  28  x  21  (x,y,A.)  voxels.  There  was  no  allowance  for 
atmospheric  absorption  over  the  path  difference  external  to  the  field  stop,  but  absorption  over  such  a  short 
path  should  not  be  significant  outside  of  the  CO2  absorption  band  from  4.2  to  4.4  pm.  Of  greater  concern 
is  the  absorption  of  the  fiber  and  the  spectral  variation  of  the  monochromator/blackbody  output.  These 
factors  will  be  corrected  in  future  calibrations  by  use  of  a  reference  detector  at  the  focal  plane  of  the  Ge 
minifying  lens. 

It  is  not  yet  established  whether  the  objective  lens  is  exchangeable,  as  in  the  visible  instrument, 
under  these  calibration  conditions.  Assuming  low  aberrations  from  2.95  to  5.05  pm,  (especially 
distortion)  in  the  objective  lenses,  it  should  be  possible  to  exchange  objectives  without  losing  calibration 
accuracy.  It  is  very  simple  to  factor  out  the  ratio  of  the  spectral  transmission  of  the  exchanged  objectives 
from  the  calibration  images. 


RESULTS 


Preliminary  results  have  been  obtained  using  two  laboratory  targets.  The  first  was  a  rear 
illuminated  filter  wheel  with  an  open  position  and  a  3.4  -  4.5  pm  bandpass  filter.  A  lighter  produced  hot 
combustion  gases  which  can  be  seen  in  the  lower  right  comer  of  the  raw  focal  plane  image  shown  below 
in  Fig.  7.  Four  dispersed  orders  can  be  clearly  seen,  leaving  much  of  the  FPA  area  unused. 


Fig.  7:  Raw  focal  plane  snapshot  image,  with  background  subtraction, 
of  a  filter  wheel  with  a  3.4  -  4.5  pm  bandpass  filter  in  the  upper  left 
position  and  no  filter  in  the  lower  right  position.  A  hot  plate  was  used 
for  backlighting.  Hot  combustion  products  from  a  butane  lighter  can  be 
seen  in  the  lower  right  comer  of  the  center  order.  Spectral  separation  of 
the  combustion  gas  from  the  heated  plate  is  evident  in  all  but  the  center 
order. 


The  four  primary  orders  and  the  center,  spectrally  compressed,  order  were  used  to  reconstract  the  object 
via  the  iterative  Multiplicative  Algebraic  Reconstmction  Technique  (MART).^  The  grayscaled  sum  of 
all  21  spectral  images  is  shown  in  Fig.  8  and  the  spectra  of  selected  resels  are  shown  in  Fig.  9  below. 


Fig.  8:  Spectral  images  of  the  filter  wheel  and  hot  gas  for  all  21  bands 
of  the  instrument  compressed  into  a  broadband  image  to  show  spatial 
resolution.  The  image  is  28  x  28  resolution  elements  (resels),  Resels  A 
and  B  have  been  marked  for  further  examination. 


Fig.  9:  Spectra  of  resels  designated  A  and  B  in  Fig.  8.  Resel  A  shows 
the  influence  of  the  3.4  ~  4.5  |im  bandpass  filter  and  resel  B  is  dominated 
by  hot  CO2  emission  from  the  butane  lighter.  The  relative  brightness  of 
the  two  resels  is  reflected  in  the  amplitudes  of  their  spectra  in  agreement 
with  Fig.  8. 


The  second  laboratory  target  was  a  cup  of  hot  coffee  spot  painted  with  a  3.39  |xm  HeNe  laser. 
The  coffee  level  in  the  cup  is  apparent  in  the  center  order  but  the  laser  spot  is  partially  masked  by  the  hot 
liquid.  In  the  four  primary  dispersed  orders,  the  laser  spot  stands  out  since  the  system’s  dispersion  places 
longer  wavelengths  farther  from  the  center  order. 


Fig.  10:  Raw  focal  plane  snapshot  image,  with  background  subtraction, 
of  a  cup  of  hot  coffee  painted  by  a  3.39  jim  laser  in  the  center  warm  area. 
Although  the  laser  spot  is  not  easy  to  discern  in  the  center  broadband 
order,  spectral  separation  of  the  spot  is  obvious  in  the  four  primary 
dispersed  orders  on  the  diagonals. 


Because  of  its  ability  to  resolve  sharp  spectral  features  such  as  emission  lines,  the  Expectation 
Maximization  (EM)  technique*  was  used  to  reconstruct  this  test  target.  Results  are  shown  in  Fig.  1 1  and 
Fig.  12  below. 
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Fig.  11:  Reconstructed  spectral  images  of  the  coffee  cup  and  3.39  pm 
laser  spot  for  all  21  bands  of  the  instrument  compressed  into  a 
broadband  image  to  show  spatial  resolution.  The  image  is  28  x  28 
resolution  elements  (resels).  Resels  A  and  B  have  been  marked  for 
further  examination.  Note  the  hot  liquid  level  is  consistent  with  Fig.  10. 
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Fig,  12:  Spectra  of  resels  A  and  B  in  Fig.  1 1  above.  Resel  A  exhibits  a 
strong  laser  signal  in  the  band  centered  at  3.4  pm  whereas  resel  B  does 
not.  Both  resels  exhibit  spectral  consistency  beyond  4  pm. 


DISCUSSION 


In  Figs.  7-12,  basic  functionality  of  the  instrument  has  been  verified.  Salient  spectral  features 
of  both  targets  have  been  successfully  resolved.  The  bandpass  filter  and  the  hot  CO2  from  combustion 
were  reconstructed  at  the  proper  spectral  wavelength  regions  in  the  MWIR  for  the  first  target.  Likewise, 
the  3.39  |i.m  laser  line  peak  appeared  in  the  spectral  reconstruction  band  centered  at  3.4|i.m  for  the 
second  target.  Both  spectrally  compressed  reconstructions  in  Fig.  8  and  Fig.  1 1  demonstrated  sufficient 
spatial  resolution  to  identify  pertinent  features  of  the  original  scene.  For  instance,  the  filter  wheel 
apertures  in  Fig.  8  and  the  coffee  cup  handle  in  Fig.  1 1  are  clearly  identifiable.  The  instrument  is 
presently  limited  by  the  first  generation  binary  diffraction  grating  which  poorly  utilizes  focal  plane  array 
real  estate  and  places  a  disproportionately  large  amount  of  light  into  the  center  order.  A  second 
generation  holographic  disperser,  which  will  produce  25  uniform  orders,  is  being  fabricated  to  rectify 
these  problems. 


CONCLUSION 

We  have  demonstrated  a  high-speed  snapshot  imaging  spectrometer  operating  in  the  MWIR. 
Relative  spectral  mapping  of  an  object  has  been  demonstrated,  although  an  absolute  radiometric 
calibration  of  the  instrument  has  not  yet  been  performed.  Future  work  includes  performing  this 
calibration  to  enable  absolute  temperature  mapping  and  absolute  spectral  characterization. 
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Complete  Imaging  Stokes  Polarimeter  for 
The  3-5  |im  Spectral  Band 


INTRODUCTION 

A  MWIR  imaging  polarimeter  has  been  under  development  in  parallel  with  the  MWIR  CTIS. 
Measurements  of  the  state  of  polarization  of  radiation  are  of  interest  in  remote  sensing  and 
reconnaissance.  Since  radiation  emitted  or  reflected  by  smooth  surfaces  tends  to  take  on  a  polarized 
component,  polarimetry  has  the  potential  of  aiding  detection  and  identification  of  manmade  objects  in 
natural  backgrounds,  even  when  those  objects  are  spectrally  camouflaged. 

The  state  of  polarization  is  often  described  by  the  Stokes  vector  s.  Its  first  components  So 
represents  the  total  radiance  (or  other  radiometric  derivative  of  power)  of  a  pixel  in  the  target  scene.  The 
latter  three  components  (s/,  S2,  and  S3)  represent  respectively  preferences  for  horizontal  vs.  vertical  linear 
polarization,  linear  polarization  at  45°  vs.  135°,  and  right  vs.  left-handed  circular  polarization.  A  number 
of  other  quantities  for  describing  the  state  of  polarization  can  be  defined  in  term  of  the  Stokes  vector.’ 


These  include  the  degree  of  polarization  (DOP),  the  degree  of  linear  polarization  (DOLP),  and  the 
orientation  angle  7]  of  the  major  axis  of  the  polarization  ellipse. 


DOP  =  (j,  +S2+  )/  Sq 

(1) 

DOLP  =  (5,  +  ^2 )/  •^0 

(2) 

77  =  4arctan(52  /5,) 

(3) 

The  MWIR  imaging  polarimeter  is  in  particular  a  complete  Stokes  polarimeter,  meaning  it  is 
capable  of  measuring  all  four  components  of  the  Stokes  vector  s.  It  features  a  rotatable  retarder  and  fixed 
linear  polarizer  in  the  fore  optics  (an  architecture  which  will  be  abbreviated  RRFP  hereafter).  Images  are 
captured  with  the  retarder  rotated  about  the  optical  axis  to  each  of  several  predetermined  angles,  and  a 
Stokes  vector  estimate  for  each  pixel  in  the  images  is  reconstructed  from  this  raw  data.  The  choice  of 
values  of  retardance  and  retarder  orientation  angles  was  made  on  the  basis  of  the  noise  sensitivity  in  the 
reconstruction  step,  as  will  be  discussed  shortly. 


SYSTEM  HARDWARE 

Figure  1 3  contains  a  representation  of  the  polarization  and  imaging  optics  in  the  polarimeter.  A 
wire  grid  polarizer  (The  Optometries  Group,  Cat.  No.  5-2625)  and  4.4  -  5.5  |i.m  bandpass  filter  are  located 
together  in  a  filter  wheel  position  inside  the  camera  dewar.  The  retarder  is  a  form-birefringent  grating* 
fabricated  in  GaAs  at  Sandia  National  Laboratories  (Fig.  14).  It  has  a  nominal  retardance  of  120°  in  the 
system’s  spectral  band  with  was  designed  to  for  minimal  wavelength  dependence  in  retardance.  The 
retarder  exhibits  diattenuation  with  an  extinction  ratio  on  the  order  of  1.5  which  is  attributable  to 
differences  in  Fresnel  reflection  due  to  the  differing  effective  indices  of  the  TE  and  TM  modes  of  the 
grating.  The  axes  of  the  diattenuation  are  coincident  with  the  fast  and  slow  axes  of  the  retarder.  As  a 
result  the  diattenuation  can  be  employed  to  accurately  determine  the  orientation  of  the  wire  grid  polarizer 
in  the  dewar  by  means  of  a  crossed-polarizer  effect. 

The  12  bit,  512  x  512  element  InSb  camera  system  described  in  the  discussion  of  MWIR  CTIS 
has  been  used  with  this  polarimeter.  It  affords  significant  improvement  over  the  data  obtained  during 
initial  system  tests  with  an  8  bit  PtSi  imager  reported  in  Ref.  9,  making  new  issues  of  system  calibration 
and  characterization  accessible  to  study. 


Fig.  13;  Schematic  showing  the  retarder  and  polarizer  of  the  RRFP  system  with 
an  objective  lens. 
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Fig.  14:  Part  a  shows  measured  and  simulated  (using  rigorous  coupled  wave 
analysis)  values  of  retardance  of  a  form-birefringent  retarder  as  a  function  of 
wavelength.  Part  b  is  a  scanning  electron  micrograph  of  such  an  element, 
showing  subwavelength  structure.  These  data  refer  to  a  test  piece  fabricated 
separately  from  the  retarder  used  in  the  polarimeter,  and  the  retardance  differs 
from  that  quoted  in  the  text.  See  Ref.  9  for  details. 


SYSTEM  PERFORMANCE 

The  retarder’s  retardance  and  orientation  angles  were  chosen  to  minimize  the  noise  sensitivity  of 
the  process  of  reconstructing  Stokes  vector  estimates  from  the  raw  data.  Sensitivity  to  noise  was 
characterized  by  several  previously  described  figures  of  merit'®  based  on  the  measurement  matrix’  for  an 
RRFP  system.  The  measurement  matrix  arises  in  a  general  linear  model  of  a  polarimeter  as  an  instrument 
which  acquires  data  by  measuring  the  power  remaining  in  a  beam  after  passing  through  a  system  of 
polarization  optics  described  by  known  Mueller  matrices.  If  a  total  of  N  power  measurements 
A:  e  {0,1,...A^-1 }  are  made,  each  through  a  system  of  Mueller  matrix  M*,  the  measurements  are  related  to 
the  Stokes  vector  s  of  the  incident  beam  in  the  absence  of  noise  via 


p  =  Ws. 


(4) 


Here  p  is  an  A^-element  column  vector  containing  the  pk  and  W  is  the  Nx4  measurement  matrix.  The  k'*' 
row  of  W  is  computed  as  a^Mk,  where  a  is  a  4-element  column  vector  describing  the  polarization- 
sensitive  responsivity  of  the  detector. 

An  estimate  s  of  the  incident  Stokes  vector  is  generated  as 

s  =  p .  (5) 

Here  a  superscripted  +  represents  the  pseudoinverse"  operation.  The  pseudoinverse  is  used  for 
generality,  since  it  can  be  applied  even  when  W  is  not  invertible.  In  the  case  of  N=4  measurements, 
however,  a  complete  polarimeter  will  by  definition  have  an  invertible  measurement  matrix.  In  order  to 
describe  a  real  system  in  which  noise  is  present,  Eq.  4  should  be  replaced  by 

p  =  Ws  +  n  (6) 

where  n  is  an  A^-element  noise  vector  whose  elements  may  be  considered  random  variables.  Substitution 
into  Eq.  5  yields  a  term 

£  =  W^n  (7) 

which  represents  the  error  due  to  noise  in  the  Stokes  vector  measurement.  A  number  of  figures  of  merit 
have  been  used  to  quantify  the  effect  of  the  system  matrix  W  in  this  term.  These  include  the  determinant 
of  W  and  the  ratio  of  its  largest  to  smallest  singular  values.  We  prefer  a  figure  of  merit  which  we  denote 
EWV  that  has  a  simple  interpretation  in  terms  of  the  noise  sensitivity  of  the  reconstruction  process. 
EWV,  which  stands  for  equally  weighted  variance,  can  be  expressed  in  several  forms. 

£W=XX(W^)%  =Tr{(W")"wn  =  |;i///;  (8) 

7=0  k=0  7=0 

Here  the  Pj  are  the  singular  values  of  W  and  R  is  its  rank.  (R=4  for  the  case  of  a  complete  polarimeter.) 
In  the  case  that  the  components  of  n  are  independently  random  with  the  same  variance  u,  EWV  is  equal  to 
the  sum  of  the  variances  of  each  of  the  components  of  s  normalized  by  u.  It  is  proportional  to  an  average 
(over  the  components  of  the  Stokes  vector)  of  the  gain  in  variance  incurred  in  the  reconstruction. 

Numerical  minimization  of  the  EWV  for  a  RRFP  system  (using  N=4  measurements  with  a 
polarization  insensitive  detector  with  unit  responsivity  and  optics  of  unit  transmittance)  yields  an 
optimum  value  of  132°  (approximately  3/8  wave)  for  the  retardance. The  value  of  EWV  at  the 
minimum  is  10.0.  Optimum  values  for  the  retarder  orientation  angle  for  each  measurement  were  also 
obtained  from  the  optimization.  The  resulting  configuration  is  also  optimal  with  respect  to  the 
measurement  matrix  determinant  and  the  ratio  of  its  largest  to  smallest  singular  values. 

The  foim-birefringent  retarder  was  designed  for  the  optimal  retardance  of  132°,  and  the  actual 
value  was  measured  at  116.8°.  Retarder  orientation  angles  for  use  in  practice  were  chosen  via  a  similar 
optimization  using  a  system  matrix  which  accounts  for  the  diattenuation  of  the  retarder  and  the  actual 
retardance.  The  resulting  angles  are  ±53.7°,  ±15.6°  (between  the  retarder’s  fast  axis  and  the  transmissive 
axis  of  the  polarizer).  The  EWV  of  the  measurement  matrix  when  these  angles  are  used  and  diattenuation 
is  included  is  increased  moderately  to  14.6. 

We  have  found  it  necessary  to  subtract  a  pedestal  image  from  each  raw  image  of  the  target  scene 
before  reconstructing  Stokes  vector  estimates.  Offsets  due  to  fixed  pattern  noise  in  the  camera  system 
and  narcissus  effects  in  the  optics  can  be  eliminated  in  doing  so.  Since  narcissus  effects  vary  with 
rotations  of  the  retarder,  a  separate  pedestal  image  is  necessary  for  each  retarder  orientation  angle.  Ideally 
the  pedestal  images  would  be  obtained  by  filling  the  field  of  view  with  a  cold  object  in  the  object  plane. 
This  is  sometimes  possible  in  a  laboratory  setting,  but  may  be  unrealistic  in  a  remote  sensing  scenario. 
An  alternative  is  to  capture  the  pedestals  with  a  cold  baffle  immediately  in  front  of  the  system’s  aperture. 
Subtraction  of  pedestals  obtained  by  this  method  does  not  correct  for  atmospheric  emission  over  the 


viewing  path,  however,  which  may  be  important  depending  on  the  spectral  band  in  use  and  the  viewing 
distance. 

The  spectral  band  of  this  instrument  was  chosen  on  the  long  wavelength  side  of  the  CO2 
absorption  feature  in  the  3-5  |a,m  atmospheric  transmission  window  in  order  to  minimize  sensitivity  to 
reflected  solar  radiation  (and  thereby  isolate  emitted  radiation).  A  LOWTRAN'^  calculation  suggests 
that,  depending  on  humidity,  atmospheric  emission  will  make  a  non-negligible  contribution  to  our 
system’s  measurements  in  the  4.4-5. 5  |im  band  over  viewing  distances  on  the  order  of  100  m.  The 
atmospheric  emission  will  effectively  be  a  superposition  of  a  uniform,  randomly  polarized  field  onto  the 
scene,  i.e.  an  addition  to  the  So  Stokes  component.  The  DOP  and  DOLP  will  be  reduced  as  a  result. 
However  the  polarization  angle  77,  which  is  independent  of  So,  will  be  unaffected  and  may  therefore  still 
be  a  useful  parameter  for  detection  and  identification  of  targets.  Nonetheless,  it  may  be  necessary  to 
reduce  the  system’s  spectral  bandwidth  or  move  the  band  to  the  short  wavelength  side  of  the  CO2 
absorption  feature  to  reduce  the  effects  of  atmospheric  emission. 


CONCLUSION 

Careful  treatment  of  such  concerns  as  atmospheric  emission  are  necessary  before  precise  Stokes 
vector  measurements  are  possible.  Fig.  shows  qualitative  results  of  a  laboratory  test  in  which  the  target 
was  a  wire  grid  polarizer  mounted  in  a  rotation  stage  in  front  of  a  wide-aperture  blackbody.  Specular 
reflections  of  the  environment  from  the  face  of  the  wire  grid  polarizer  were  blocked  with  a  cryogenically 
cooled  baffle  during  image  capture.  The  face  of  the  polarizer,  which  was  oriented  to  transmit  horizontally 
polarized  radiation,  is  prominent  in  the  Si  and  DOLP  images. 


Si 

Fig.  15;  So,  Si,  and  DOLP  images  from  the  experiment  described  in  the 
text.  The  5/  and  DOLP  images  are  displayed  in  reverse  contrast. 
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